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Certain two-component reaction-diffusion systems on a finite interval are known to possess mesa (box-like) steady- 
state patterns in the singularly perturbed limit of small diffusivity for one of the two solution components. As the 
diffusivity D of the second component is decreased below some critical value D c , with D c = O(l), the existence of 
a steady-state mesa pattern is lost, triggering the onset of a mesa self-replication event that ultimately leads to the 
creation of additional mesas. The initiation of this phenomena is studied in detail for a particular scaling limit of the 
Brusselator model. Near the existence threshold D c of a single steady-state mesa, it is shown that an internal layer 
forms in the center of the mesa. The structure of the solution within this internal layer is shown to be governed by 
a certain core problem, comprised of a single non-autonomous second-order ODE. By analyzing this core problem 
using rigorous and formal asymptotic methods, and by using the Singular Limit Eigenvalue Problem (SLEP) method 
to asymptotically calculate small eigenvalues, an analytical verification of the conditions of Nishiura and Ueyema 
[Physica D, 130, No. 1, (1999), pp. 73-104], believed to be responsible for self-replication, is given. These conditions 
include: (1) The existence of a saddle-node threshold at which the steady-state mesa pattern disappears; (2) the 
dimple-shaped eigenfunction at the threshold, believed to be responsible for the initiation of the replication process; 
and (3) the stability of the mesa pattern above the existence threshold. Finally, we show that the core problem is 
universal in the sense that it pertains to a class of reaction-diffusion systems, including the Gierer-Meinhardt model 



with saturation, where mesa self-replication also occurs. 

X 



1 Introduction 

In [28] Pearson used numerical simulations to show that the two-component Gray-Scott reaction-diffusion model 
in the singularly perturbed limit can exhibit many intricate types of spatially localized patterns. Many of these 
numerically computed patterns for this model have been observed qualitatively in certain chemical experiments 
(cf. [12] . |13| ). An important new phenomenon that was discovered in [28] . [12] , and [13] , is the occurrence 
of self-replication behavior of pulse and spot patterns. In recent years, many theoretical and numerical studies 
have been made in both one and two spatial dimensions to analyze self-replication behavior for the Gray-Scott 
model in different parameter regimes (cf. [33], [32], [26], [27], [34], [20], [2], pQ, [16]). In addition to the Gray- 
Scott model, many other reaction-diffusion systems have been found to exhibit self-replication behavior. These 
include the ferrocyanide-iodide-sulfite system (cf. [13] ). the Belousov-Zhabotinsky reaction (cf. [19j ). the Gierer- 
Meinhardt model (cf. [H], [4], [17]), and the Bonhoffer van-der-Pol-type system (cf. [§], [9]). 

Despite a large number of studies on the subject, the detailed mechanisms responsible for self-replication are 
still not clear. In an effort to classify reaction-diffusion systems that can exhibit pulse self-replication, Nishiura 
and Ueyema [26] (see also [5]) proposed a set of necessary conditions for this phenomenon to occur. Roughly 
stated, these conditions are the following: 
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(a) (b) 

Figure 1. (a) Numerical simulation showing mesa-splitting in the Brusselator model (|1.2| . The fixed parameters 
are /3q = 1.5, e = 0.01, r = 0.7, x E [0, 2]. The parameter D is slowly decreased in time according to the formula 
D = (1 + 5 x 10~ 6 t)~ 2 . Parameters (3 and a are determined through (|1.3|) . i.e. a = e 2 /D, (3 = a/3o- The scale on 

the vertical axis is if = , where £> c = 0.88. Splitting events occur for K w 1 and K k, 2. (b) Snapshots of 
the profile of m during a splitting event. The time between two successive snapshots is 1000 time units. 



(1) The disappearance of the if-spike steady-state solution due to a saddle-node (or fold point) bifurcation 
that occurs when a control parameter is decreased below a certain threshold value. 

(2) The existence of a dimple eigenfunction at the existence threshold, which is believed to be responsible for 
the initiation of the pulse-splitting process. By definition, a dimple eigenfunction is an even eigenfunction 
$(y) associated with a zero eigenvalue, that decays as \y\ — > oo and that has precisely one positive zero. 

(3) Stability of the steady-state solution above the threshold value for existence. 

(4) The alignment of the existence thresholds, so that the disappearance of K pulses, with K — 1,2,3,..., 
occurs at asymptotically the same value of the control parameter. 

For the Gray-Scott model in the weak interaction parameter regime where the ratio of the diffusivities is O(l), 
Nishiura and Ueyema in [26] verified these conditions numerically for a given fixed diffusivity ratio. Alternatively, 
for the Gray-Scott model in the semi-strong regime, where the ratio of the diffusivities is asymptotically large, it 
was shown in |20[ and in equation (2.9) of [2] that the following core problem determines the spatial profile of a 
pulse in the self-replicating parameter regime: 

V'-V + UV 2 =0, U"-UV 2 = 0; U' (0) = V' (0) = , V -> , U' -» A as y -> oo . (1.1) 

By using a combination of asymptotic and numerical methods, and by coupling (|1.1| to an appropriate outer 
solution away from a localized pulse, conditions (l)-(4) of Nishiura and Ueyema were verified in |16) . In [3j a 
detailed study of the intricate bifurcation structure of (|1.1|) was given. 

In this paper we study self-replication of mesa patterns. A single mesa solution is a spatial pattern that consists 
of two back-to-back transition layers. An example of such a steady-state pattern is shown in Fig. [2] below. Our 
goal is to analytically verify whether the conditions (l)-(4) of Nishiura and Ueyama [26], originally formulated 
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for analyzing pulse self-replication behavior, also hold for mesa self-replication. In addition, we seek to derive and 
study a certain core problem, analogous to that pertains to self-replicating mesa patterns. 

For concreteness, we concentrate on the Brusselator model. This model was introduced in [30] . and is based on 
the following hypothetical chemical reaction: 

A -> X , C + X Y + F , 2X + Y ->3X, X -> E. 

The autocatalytic step 2X + Y — > 3X introduces a cubic non-linearity in the rate equations. Since the 1970's, 
various weakly-nonlinear Turing patterns in the Brusselator have been studied both numerically and analytically 
in one, two, and three dimensions. These include spots, stripes, labyrinths and hexagonal patterns (cf. [6], |21j . 
[29] . [35] . [36]), and oscillatory instabilities and spatio-temporal chaos (cf. |14j . |37| ). 
After a suitable rescaling, we write the one-dimensional Brusselator model on a domain of length 2L as 

u t — e 2 u xx - u + a + u 2 v , rv t = e 2 v xx + (1 - 0) u - u 2 v ; u x (±L, t) = V X (±L, t) = . (1.2) 

In this paper we make the following assumptions on the parameters: 

e 2 

£<1; a<l; < 1 ; D = — = O(l) ; #, = - = 0(1) , with O > 1 ; r = 0. (1.3) 

a a 

The full numerical results in Fig. Q] illustrate the mesa self-replication behavior for (|1.2|) . To trigger mesa self- 
replication events we started with a single mesa as initial condition and slowly decreased D in time (see the figure 
caption for the parameter values). At the critical value D\ ~ 0.8, a mesa splits into two mesas, which then repel 
and move away from each other. The splitting process is repeated when D is decreased below D2 ~ 0.2. 

In §2 we calculate a threshold value D c of D for the existence of a single-mesa steady-state solution for the 
Brusselator (|1.2j) in the limit e — > 0, and under the assumptions (|1.3|) on the parameter values. The result, 
summarized in Proposition Q] of §2, shows the existence of a value D c such that a if-mesa steady-state solution 
exists if and only if D > D c /K 2 . Analytical upper and lower bounds for D c are also derived. Similar thresholds for 
the existence of steady-state mesa patterns were derived for other reaction-diffusion systems in [llj using more 
heuristic means. Our analysis is based on a systematic use of the method of matched asymptotic expansions. 

For a single-mesa steady-state solution, we show in §2.1 that an internal layer of width 0(e 2 / 3 ) forms in the 
center of the mesa when D is asymptotically close to the threshold value D c . This internal layer is illustrated 
below in Fig. [4j By analyzing this internal layer region using matched asymptotic analysis, we show that the 
solution u is determined locally in terms of the solution U(y) to a single non-homogeneous ODE of the form 

U" = U 2 -A-y 2 ; U'(0) = 0, U'(y)^l as y -> 00 . (1.4) 

Here A is related to the parameter values in (|1.2j) . We refer to (|1.4|) as the core problem for the onset of self- 
replication. Unlike (|1.1|) for self-replicating pulses in the Gray-Scott model, the problem (|1.4j) is not coupled and, 
consequently, is easier to study analytically than (|1.1|) . The proof of conditions (1) and (2) of Nishiura and Ueyama 
is then reduced to a careful study of (II. 4|) . More specifically, by using rigorous techniques we prove analytically the 
existence of a saddle-node bifurcation for (II. 4p and we analyze the solution behavior on the bifurcation diagram. 
The result is summarized below in Theorem[2] In §2.2 we use some rigorous properties of the core problem, together 
with a formal matched asymptotic analysis, to asymptotically construct a dimple eigenfunction corresponding 
to the zero eigenvalue at the saddle-node bifurcation value. This construction, summarized in Proposition [31 
establishes condition (2) of Nishiura and Ueyama. 

In §2.3 we show that the core problem (|1.4|) is universal in the sense that it can be readily derived for other 
reaction-diffusion systems where mesa self-replication occurs. The universal nature of (II. 4|) is illustrated for some 
specific systems, including the Gierer-Meinhardt model with saturation (cf. |18j ). For this specific model, mesa- 
splitting was computed numerically in Figure 28 of [17] . Although the phenomena of mesa self-replication is 
qualitatively described in Chapter 11 of [11] . the core problem (II. 4|) governing the onset of mesa self-replication 
and its analysis has not, to our knowledge, appeared in the literature. 

Since the saddle- node existence value for a if- mesa steady-state solution is D = D c /K 2 , the condition (4) of 
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Xi- Xi + x 2 - x 2 + 

-L -I I L 

(a) (b) 

Figure 2. (a) A single mesa steady-state on [—L,L]. The parameter values are L = 1, /3o = 1.5, D — 1.29, e = 
0.005. (b) A two-mesa steady state [K = 2, L = 1/4, x £ [0,1]). The parameter values are D = 0.068, e = 
0.00125, /? = 1-5. 

Nishiura and Ueyama, regarding an asymptotically close alignment of saddle-node bifurcation values, does not 
hold in a strict sense. However, this condition is satisfied in the same approximate sense as in the study of self- 
replicating pulses for the Gray-Scott model in the semi-strong interaction regime (see Table 3 of pQ and equation 
(1.2) of [16]). 

In $3] we study the stability of if- mesa steady-state solutions when D > D c /K 2 . We show that such a pattern 
is stable when r = 0, and moreover all asymptotically small eigenvalues are purely real. This proves condition (3) 
of Nishiura and Ueyema. Our analysis is similar to the SLEP method, originally developed by Nishiura et. al. in 
[22] and [23j , and that has been used successfully to prove the stability of mesa- type patterns in reaction-diffusion 
systems and in related contexts (cf. [24 , [25 ). In our analysis a formal matched asymptotic analysis is used to 
derive a reduced problem that capture the asymptotically small eigenvalues of the linearization. This reduced 
system is then studied rigorously using several tools, including the maximum principle and matrix theory. In this 
way we prove that the small eigenvalues are purely real and negative when r = 0. 

Finally, in Sj4] we relate our results regarding mesa self-replication for the Brusselator model with previous 
results concerning the coarsening phenomenon of mesa patterns that occurs when D is sufficiently large (cf. |15| ). 
In addition, we propose some open problems. 

2 The Steady State Mesa and the Universal Core Problem 

In this section we study the steady-state problem for (|1.2|l . and we prove analytically the first two conditions of 
Nishiura and Ueyema. We will analyze an even symmetric solution of the type shown on Fig. [21(a), consisting of 
a single mesa on a domain [-L, L] with interfaces at x = ±L The K-mesa solution on a domain of length 2KL 
can then be constructed by reflecting and gluing together K such solutions. 
We first reformulate (|1.2|) to emphasize the slow-fast structure. We define w by 



W = V + u , 
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so that (|1.2|) becomes 

Ut — e 2 u" — u + a + u 2 (w — u) , r (wt — u t ) + Ut — e 2 w" + a — 8 . (2-1) 

Here the primes indicate derivatives with respect to x. We then introduce 8q = 0(1) and D = 0(1) defined by 

8 e 2 
8 = ^-, D =— . 
a a 

Then, (|2.ip becomes 

t (1 — rl 

u t = eV - u + a + u 2 (w - u) , —w t +- -u t = Dw" + 1 - /3 . (2.2) 

a a 

The corresponding steady-state problem is 

s 2 u" -u + u 2 (w~u) + a = 0, Dw" + l-8 u = 0, (2.3) 

where a = 0(e 2 ). Since e 2 -C D from (| 1 . 3[) . it follows that w is the slow variable and u is the fast variable. Upon 
integrating (|2.3|) for w and using the Neumann boundary condition for w at x = L, and the symmetry condition 
m/(0) = 0, we obtain the integral constraint 



L = B udx. (2.4) 
Jo 

Near the interface at x = I we introduce the inner expansion 

u = U (y)+eU 1 (y) + w = W + eWi (y) + . . . , y = e" 1 ^ - /) . (2.5) 
Upon substituting this expansion into (|2.3|) . we obtain the leading-order problem 

U5'-/(l/o J Wb) = J < = 0, (2.6) 

where f{u, w) is defined by 

f(u,w) = u- u 2 (w - u) . (2.7) 

At next order, we obtain 

CUt = U'( - f u (U Q , W )Ui = f w {U ,W )W x , W'{ = . (2.8) 

From ()2.6j) we get that Wo is a constant to be determined. To ensure that there exists a heteroclinic connection 
for U~o we require that / satisfy the Maxwell line condition, which states that the area between the first two 

roots of / is the negative of the area between its last two roots of /. Since / is a cubic, this is equivalent to 
simultaneously solving / = and /" = for Wq- In this way, we obtain 



1 ± tanh ( - 



(2.9) 



For the mesa solution as shown in Fig. [2ja) , we must take the minus sign in (|2.9[) above. 

To determine the interface location I, we now study the outer problem away from the interface at x = I. Since 
a = 0(e 2 ), we obtain to leading order from (|2.3[) that 

u + u 2 (w — u) = . 

This yields either u = or 

w~h(u) = -+u. (2.10) 

u 

Moreover, we have Uq — * as y — > oo and Uq — * V2 as y — > — oo. Therefore, by matching to Uq and Wo, and by 
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using the symmetry condition at x — 0, we obtain the following outer problem in the mesa region < x < I: 

w = h{u); Dw" = g{u) = (3 u- 1, < x < I. (2.11a) 

u(l) = V2, w{l) = ^=, u' (0) = w' (0) = . (2.116) 
v2 

The leading-order outer problem on I < x < L \s u = and Dw" = — 1. 

The solution to the second-order inner problem (|2.8|) for W\ is W\ = W\\y + Wvi, where W\\ and W12 are 
constants to be determined. Since CUq = 0, the solvability condition for (|2.8p yields 

/oo />oo 
U^f w (U , W )W 1 dy = - U' a U% (Wny + W 12 ) dy . 
-00 J —oo 

This yields one relation between W\\ and W\i. The second relation is obtained by matching W to the the outer 
solution w. This yields Wn = w'(l ). In this way, we obtain 

W[ = W 11 =w'{l), W 12 = -^±J y(U*)' dy. (2.12) 

We now solve the outer problem (|2. 1 1[) in terms of uq = u(0). We first define F (u; uq) by 

F(u;u )= f g(s)h'(s)ds. (2.13) 

Ju 

By multiplying (|2.11 al) for w by w' we get 

D—=F(u;u ), w' = J — y/F{u;u ). 



In the outer region on / < x < L, we have u — 0. Therefore, by integrating u>" from a; = to x = L, we obtain 
/ ' #(w) dx + j. L (—l) c£x = 0. This yields, 

1 

g(u)dx = L-l. (2.14) 



The left-hand side of (|2.14p is evaluated by integrating w" from x = to x — I to get J* g (u) dx — Dw' (I) 



(2D)F(y/2;u ). In addition, by using w' = h'(u)u' , we obtain 



du / 2F (u; uq) , , 



dx V D 

We then integrate (|2.15p with u(0) = uq and u(l) = In this way, we obtain 

,s/2 



[h'iuT 1 . (2.15) 



Upon integrating the second expression in (|2.16p by parts we get 



a/2F («;u ) .9(V2) J„ [g(uf 

By combining this relation with (|2.16p . and by calculating g(\/2), we obtain the following expression relating uq 
to the overall length of the domain: 

xM^vW^o) h^7rV)+ r frL^F{u-,u )du = ±=. (2.i7) 

We note that the function ft, (u) has a minimum at uq = 1, and that = — g (uq) /i' (uq) < for > 1 from 



Self- Replication of Mesa Patterns in Reaction- Diffusion Systems 7 

(|2.13|) . Hence, \ ( u o) is a decreasing function of uq for uo > 1. Upon defining D c by D c — L 2 j\ (l) 2 , we conclude 
that there exists a value of Uo, and consequently an outer solution on < x < I exists, if and only if D > D c . 
Therefore, we obtain a threshold value D c for the existence of a single mesa solution on [— L, L\. 

Next, we obtain explicit bounds on D c . The function w = h (u) satisfies h (1) = 2 and h (V2) = 3/v2- Therefore, 
for x € [0, 1] we have that u and w are both increasing functions for x € [0, 1] with u (I) = y/2, w (I) = 3/\/2, and 

Ku<V2, 2<w<3/v / 2. (2.18) 

Using f2~T5|) and (|2.11 al) for w we estimate w" = D" 1 (/3 u - 1) > D^iPo - 1)- This yields that 

> w(0) + (/3 °~ 1) ^ 2 , which implies > 2 + ^--^-l 2 . (2.19) 

v y ~ w 2£> V2 ~ 2D v y 

Then, we use (|2.4p with u = for Z < a; < L to get 

L = (3 Q [ udx < /3 lV2. 
Combining this inequality with (|2.19[) , we get 

I 2 > — , -= > 2 1 



2/3 2 ' V2" 4/3 2 £>' 
By solving this second relation for D we conclude that a necessary condition for the existence of a solution to 
(|2.11|) when e < 1 is that D > D, where 



D= ^°- 1) I L 2 . (2.20) 
- 4/3 2 3/V2-2 



This relation provides a lower bound for D c . To find an upper bound we calculate w" < (f3 V2- l) /D. This 
inequality can be integrated, and with w(l) = V3/2, we get 



V2 2D 

Then, from (|2.4[) . we obtain 



l<±-. (2.21) 

Po 

The condition w (0) > 2 is satisfied provided that — jp-^z (pov^ - l) > 2. Hence, a sufficient condition for 
the existence of a solution to (|2.11[) when e <C 1 is that D > D, where 

D=^ 1 r 

2fi 2 3/V2-2 

It follows that D < D c < D. We summarize the results of this analysis of a single mesa solution as follows: 

Proposition 1 (Nishiura-Ueyema's Condition 1: The Steady State and its Disappearance) Consider the steady 
state solution to the Brusselator i2.S\) with (3q > 1 in the limit e — > 0. We define F(u;uo) and x( u o) by 

F(u;u Q )= (A)S-1) U- — )ds, x(uo) = \J2F(V2;uo) \ — - +/ — du . 

Ju \ s 2 J v \PoV2-lJ J U0 (Pou-iy 

In terms o/x(l), we define the threshold D c by 

D C = L 2 / X (lf. 

Suppose that D > D c . Then, there exists a uq G (l, V2) and I £ (0, L) given implicitly by 



xM = -^=, l = L-VDy/2F(V2;uo) 
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Figure 3. The graphs of D c , D and D_ versus (3$ 

such that there exists a symmetric mesa solution on the interval [~L, L] with interfaces at ±1 and with u (0) = uq. 
In the region x € (0, 1), w and u are given implicitly by 



w=- + u, Dw" — g (u) ee (3 Q u - 1 , < x < I ; u(Q)=u , u(l) = V%- 
u 

In the region x S {I, L), the leading-order outer solutions for u and w are 

« = 0, W ~-^(x-L) 2 + ^(2/V3-l) 2 + 2/V3. 



(2.22) 



Moreover, we have D < D c < D where 

(Po - 1) 



D 



4/3 2 3/V2-2 ' 



Tr = V2/3 -l 1 l2 
2PI 3/V2 - 2 



The graphs of D c , D, and D_, versus @o are shown in Fig. A single-mesa solution does not exist if D < D c . By 
reflections and translations, a single-mesa solution can be extended to a K-mesa solution on an interval of length 
2KL. 



2.1 Core problem 

In this section we verify Nishiura and Ueyema's second condition by analyzing a limiting differential equation 
that is valid in the vicinity of the critical threshold D = D c . When D is decreased slightly below D c (at which 
point u(0) ~ l,u>(0) ~ 2) the single-mesa solution ceases to exist. To study the solution near this fold point, we 
fix u(0) = uq and consider D = D(uq). From numerical computations of the steady state solution as shown in 
Fig. 21 an internal layer forms near the origin when Uq is decreased below u(0) = 1. 

To study the initial formation of this internal layer near the origin, we expand (|2.3[) near D = D c as 



u = 1 + 5u\ 



w = 2 + 8 2 wi 



D = D r 
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Figure 4. Formation of a boundary layer near the center of a mesa. The steady state and D are solved simulta- 
neously, while u(0) is fixed at one of the values 1.4, 1.3, 1.2, 1.0, 0.9, 0.4, 0.1. The corresponding values for D are as 
follows: u(0) = 1.4,13 = 17.109; u(0) = 1.3, D = 2.257; u(0) = 1.2, D = 1.290; u(0) = 1.1, L> = 0.96; u(0) = 
1.0,13 = 0.846; u(0) = .9, D = 0.863; u(Q) = .8, D = 0.938; tt(0) = .7, -D = 1.068; u(0) = .6, Z> = 1.27; u(0) = 
.5,13 = 1.624; u(0) = .4, Z3 = 2.244; u(0) = .3, 13 = 3.525; «(Q) = .2, D = 6.982; u(0) = .1, -D = 24.34. 

The nonlinear term in (|2.3p becomes f(u, w) = <S 2 (u 2 — w±) H . From (|2.3|) we then obtain 

r 2 

— < = u 2 - tui , i3 c (5 2 wi' = /3 - 1 • (2.23) 
o 

We introduce the inner-layer variable z by z = x/8 with i5< 1. Then, (I2.23P for w\ becomes 

A)-i 

The solution is 

ioi = A + Bz 2 , 6 = ^— - > , ,4 = wi (0) . (2.24) 

Then, (|2~^3]l for m is 

=u 2 - (.4 + Sz 2 ) . (2.25) 

This suggests the internal-layer scaling S — e 2 ^ 3 so that u\ zz — u\ — LA+ Bz 2 ) . The boundary conditions for u\ 
are iti, (0) = and u x -> zVB as z — ► oo. Finally, we introduce £/ and j/ as 

ui=S 1/3 C7, z = ZT 1/6 y. 

This yields the following core problem for £/(?/) on < y < oo: 

f/" = C/ 2 -A-y 2 ; C/'(0) = 0, J7' ~ 1 as y -> oo . (2.26) 

Here A is related to A and B by 

A = 2/3 . 
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Figure 5. Plot of A versus s = U(0) showing the fold point for the core problem (|2.26|) . The inserts show the 
solution U versus y at the parameter values as indicated. The middle insert shows U at the fold point. The dotted 
lines are the limiting approximations of A versus C/ ± (0) in (|2.28p . See Theorem [2] for more details. 

In terms of the original variables, we have that 

u(x) - 1 ~ e 2/3 B 1/3 U(y) , w(x) - 2 ~ s 4/3 B 2/3 (A + y 2 ) y = x/{e 2/3 B- 1/6 ) , (2.27 a) 

ui(0) = B 1/3 U(0) = e- 2/3 [«(0) - 1] , u>i(0) = B 2/3 A = e" 4 / 3 [w(0) - 2] . (2.27 6) 
Here B is defined in (12.24)) . The following main result pertains to the solution behavior of (|2.26[) . 

Theorem 2 Suppose that i> 1. Then, the core problem h2.26\) admits exactly two solutions U (y) with U' > 
for y > 0. They have the following uniform expansions: 



U+ - ^/A + y 2 , J7+(0)~\/I, (2.28 a) 

U~ - ^A + y 2 ^1 - 3 sech 2 (^^^j , U~ (0) ~ -2VA . (2.28 6) 

These two solutions are connected. For any such solution, let s = U(0) = Uq and consider the solution branch 
A = A (s). Then, A(s) has a unique (minimum) critical point at s — s c , A = A c . Moreover, define <&(?/) by 

dU 



as 



(2.29) 



Then, $ > for ally > and $ — > as y -— > 0. Numerically, we calculate that A c ps —1.46638 and s c w —0.61512. 
The graph of A (s) is shown in Fig. [3 



Proof. The proof consists of four steps. In Step 1 we use formal asymptotics to show that when A 3> 1, there 
are exactly two possible solutions with U' > for y > 0, as given by (|2.28|) . In Step 2 we rigorously show that 
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there are no solutions when —A is large enough. In Step 3 we show that the solution branch with U' > 0, y > 
cannot connect with any branch for which U' < at some y > 0. It then follows that the two branches U + and 
U~ must connect to each other. In Step 4 we show that $ is positive and that the fold point s c is unique. 

Step 1:. We first consider the case A ^> 1. After rescaling U — VAv, y = at for some a to be determined, 
(|2.26j) becomes 

1 2 « " 2 2 
=«« = « - 1 r y . 



If we choose a = \f~A then the leading-order equation for v becomes v 2 — 1 — t 2 ~ 0. This yields w ~ Vl + £ 2 - The 
other possible choice is a = A -1 / 4 , which yields the leading-order equation 

v tt = v 2 - 1 , 

with v' (0) = 0, v (£) ~ 1 for large t. This ODE admits exactly two monotone solutions satisfying v (t) — > 1 as 
t — > oo. These solutions are given by 

oft 

v = 1 and u = 1 — 3 sech 



V2 

which correspond to the inner expansion of C/ + and C/~, respectively. Matching the inner and outer expansion 
into a uniform solution yields (|2.28 a|) and (|2.28 

Step 2: Next we show the non-existence of a solution to the core problem when —A is positive and sufficiently 
large. To show this, we rescale 

u = V^Av, y=(-Ay 1/4 t. 

From (|2.26|) , we obtain 

v"=v 2 + l-et 2 , e=(-Ay 3/2 . (2.30) 

We will show that no solution to (|2.30[) exists when e > is small enough. First we choose any a € (0, 1) and 
define T by 



T = 

Then, for < t < T we have 

v" > v 2 + a ; v'(Q)=Q, v (0) = v Q . 

In particular, v' > for all £ € (0, T). First, we suppose that w (0) = Vo < 0. Then, under this assumption, wc 
derive 

w' 2 1 /l \ 1 

-y > g« 3 + a« ~ ( + a^o J > + <M ■ 

The first step is to show that when e is sufficiently large, v (£) crosses zero at some value t — t\. There two 
subcases to consider. For the first subcase, suppose that vq < — 1. Then 

v' 2 1 , 1 , 
2 - 3 3 ' 



so that 



/, / ; dv <ckr 1/2 <c 

1,3 _ 1,3 



Therefore, by choosing e small enough so that T > C, we have £i € [0,T]. For the subcase vo > —1, we have 
v" > a, v > |£ 2 + vo so that t\ < y|. Therefore, £i S [0, T] by choosing e small enough so that T > y|. 
The second step is to show that w blows up for some T\, 6 (0, T), provided that T is large enough. This would 
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yield a contradiction. Indeed we have > ^v 3 + av for v > so that 

f°° dv 
T b <h+ h, h = I , < oo. 

J o V%\J\v 3 + av 

Therefore a contradiction is attained by choosing e so small that T > I2 +t\. 

Finally, if vq > let 77 = vq + Bt 2 . Then, for large enough B we have 77" — rj 2 — 1 + et 2 < 0, so that by a 
comparison principle, v > vq + Bt 2 for all t > 0. But this is impossible since we must have v — > \fet for large 
values of t. This shows that no solution to the core problem (|2.26p can exist if —A is sufficiently large. 

Step 3: We now show that the solution branch with U' > for y > can never connect to a non- monotone 
solution branch. We argue by contradiction. Suppose not. Then consider the first parameter value A for which a 
connection occurs. For such a value of A, there must be a point yo € [0, 00) such that U' (yo) = with U' > 
for any other y. Suppose first that yo > 0. Then we have U'" = 2UU' — 2y so that V" (yo) = —2yo < 0. But this 
contradicts the assumption that U' (yo) = is a minimum of U'. If on the other hand yo = 0, then we consider three 
cases. First if U" (0) = 0, then from a Taylor expansion we obtain U (0) = VJ; U' (0) = U" (0) = U'" (0) = 0, 
(0) = —2. This expansion shows that U is decreasing to the right of the origin, which contradicts the 
assumption that U' > for all y ^ yo- Similarly, if U"(Q) < then again U is decreasing to the right of the origin, 
which yields a similar contradiction. Finally, U"(0) cannot be positive when yo = 0, since we assumed that A is 
the connection point. 

Step 4: Define $(y) by 

dU 

* = -, U(0) = s. 
At the fold point s = s c where A' (s) = 0, we obtain upon differentiating (|2.26[1 that 

$" = 2J7$, $(0) = 1. (2.31) 

To show that $ is positive at the fold point, we define x(y) by 

$ 

We readily derive that 

X "U' ~2y X + 2U" X ' = 0. (2.32) 
Since $ (0) — 1, and U' > for y > 0, we obtain that x is positive near the origin. In addition, for large y, (|2.32p 
reduces to x" ~ ^VXi which implies x ~ > as y — > 00. It follows by the maximum principle that x > 0. This 
shows the positivity of <I>. Finally, we establish that the fold point A' (s) — is unique. Assuming that A'(s) = 0., 
we differentiate (|2.26j) twice with respect to s to obtain 

= 2U<$> S + 2$ 2 - A" (s) . 

By multiplying both sides of this expression by <&, and integrating the resulting expression by parts, we obtain 

A»(s)=2 ^f* y . 
Jo 

However, since $ is positive then A" (s) > whenever A' (s) — 0. This implies that the fold point is unique. ■ 



2.2 The Dimple Eigenfunction 

Next, we study the qualitative properties of the eigenfunction pair associated with linearizing \2.2\ around the 
steady-state solution at the fold point where D = D c . We label the steady-state solution at the fold point D = D c 
by u c (x) and w c (x). From (|2.27 6|) and Theorem [21 we obtain at D = D c that 

u c0 = tt c (0) - 1 + e 2/3 B 1/3 U(0) , 17(0) = s c = -0.61512 . 
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Figure 6. (a) The plot of iti(0) versus itfi(O) computed from (I2.2T|) . The solid and dashed lines correspond 
to numerical computations of the full steady-state system (|2 . 3[) and the core problem (|2.26p , respectively, (b) 
The dimple eigenfunction at the fold point with parameter values e = 0.005, (3q = 1-5 and D is determined 
numerically as a function of ui(0) 



We linearize (I2.2|) around u c and w c by setting 



u (x, t) = u c (x) + e <p(x) , w (x, t) = w c (x) + e ip(x) . 
This leads to the eigenvalue problem 

\<j) = e 2 4> xx - f u (u c , w c ) 4>- f w (ti C! w c )ip, — [4> + t (ip - 4>)] = Dip xx - (3 (p , 



with ip x = 
Let mq 



4> x = at x = 
u(0) and D 



±L and D = D c . 

= D(uq). Then, if we define <f> and ip by 



d 



du 



u (x) 



d 



(2.33) 



(2.34) 



ip = (^) , 

it follows from ()2.33j) and D' (u c q) = that (|2.34p is an eigenfunction pair corresponding to A = 0. 

We now construct an asymptotic approximation to this eigenpair 0, ip of (|2.33j) corresponding to A = 0. In 
particular, we show that <f> is an even function that has a dimple shape when D — D c . This is shown below to 
be a consequence of the positivity of the function <i> in (12.31|) . together with the integral constraint <j)dx = 0, 
which is readily obtained from (|2.33p . We normalize this eigenfunction by imposing that </>(0) = 1. For e < 1 and 
D = D c , our analysis below shows that the asymptotic structure of <f> has four distinct regions: an inner region 
of width 0(e 2 / 3 ) near x = where cf> — O(l); an outer region on x G (0,1) where <fi = 0(e 2 / 3 ); an inner region 
of width O(e) near x = I where cf> — 0(e~ 1 ^ 3 ); and an outer region on x € [l,L] where </> = 0(e). The first three 
regions give asymptotically comparable contributions of order 0(e 2 / 3 ) to the integral constraint J Q (f>dx = 0, 
whereas the contribution from the fourth region can be neglected. For a particular set of parameter values the 
resulting dimple-shape of the eigenfunction <p at D — D c is shown in Fig. Efb) . 

We now give the details of the asymptotic construction of <fi. We begin with the internal layer region of width 
0(e 2//3 ) near x — 0. In this region, we use (|2.7p and (|2.27 a[) to calculate 



/„ - 1 + Zu 2 c - 2w c u c ~ 2e 2 / 3 B 1 / 3 U c ; f u 



Here U c (y) is the solution to the core problem (|2.26[) at the fold point location A'(s c ) = where D 
(|2.35|) in (|2 . 33[) with A = 0, wc obtain the following leading-order system on < y < oo: 



(2.35) 
Using 



£ 2 4>xx 



2e 2/3 B 1/3 U r 



■ + i> = ; D c ^ xx -M = 0, (2.36) 

with normalization condition (f)(0) = 1 and with <fi x = ip x = at x = 0. We then introduce the inner variables 
y = x/(e 2/3 B- 1/6 ), 4> = and ip = Then, (|2~36| becomes 

$"-2C/ c <I> + e" 2/3 6" 1/3 * = 0; D c ^" = /3 e 4/3 S" 1/3 <I> . (2.37) 
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This shows that 4" = 0(e 4 / 3 ), and hence e -2 / 3 '!' = 0(e 2 / 3 ) <c 1. Therefore, to leading-order, we obtain that «3> 
satisfies (|2.3ip of Theorem [2] at s = s c , and that 

* ~ ^- £ 4 /3 B -i/3^ . = ^(0) = 0. (2.38) 

Therefore, as y — > oo, we have \&o ~ (J °° <&(?/) dy) . Writing the far-field expansion of 'J in terms of the outer 
x-variable, we obtain the far-field matching condition 

*~ ^eWB- l l & (J™${y)dy\x. (2.39) 

Next, we analyze the inner region of width 0(e) near the transition layer at x = I. We introduce the inner 
variables yi — (x — l)/e, <p — $>i(y) and ip = *Bi(y). From (|2.33|) with A = 0, we obtain on — oo < yi < oo that 

(C7b, Wo) $j-/™(E7b, Wo) #i = 0, *{' = 0. (2.40) 

Here Uq and Wq, given in (|2.9p , satisfy the leading-order steady-state inner problem (|2.6p . Upon comparing (12. 6p 
and (|2.40l) . we conclude that is proportional to Uq and that 'J'; = 0. By using (12. 9|) for C/q, we get 

^-^csech 2 (^r) . ^~*; = 0. (2.41) 

Here c is an unknown constant, possibly depending on e, that is found below by the global constraint J Q L (pdx = 0. 

In the outer region x E (l,L], where u c = 0, we obtain the leading-order solution <j> = ip = 0. A higher-order 
construction, which we omit, shows that ip = cf> = 0(e) in this near-boundary region. In contrast, in the outer mesa 
plateau region x € (0, 1), we set A = in (|2.33|) to obtain <p = — [f w (u c , w c ) / f u (u c , w c )] ip and Dip xx — flo<p = 0. 
Then, by using the solution u c and w c to (|2.11[) at D = D Cl we obtain f w /f u — — u 2 /(it 2 — 1). The boundary 
conditions for ip as x — > and at x = I are obtained by matching to (|2.39[) and (|2.41|) . respectively. In this way, 
we obtain the following formulation of the leading-order outer problem for ip and <p on x G (0, 1): 

^0x^-4^-^0=0, ar€(0,I); D c ip 0x -> ft as x -> 0+ , Vo(0=0. (2.42a) 
u 2 - 1 

In terms of ipo , we have 

^ ~ e^B- 1 " (7°° t>d y ^Jip , p~ eVijL-B-^ QT $ ^ ^ . (2 . 42 6) 

By the maximum principle the solution ipo to (|2.42 ap . which depends only on /3q, satisfies ipo > on x € (0,1). 
Therefore, since <f> > from Theorem[5J we conclude that <p — 0(e 2 / 3 ) with <p > on x € (0, I) and <p(l) = 0. 

Finally, we use the global condition f Q <pdx = to calculate the constant c in (|2.41[) . Upon using <p ~ $ for 
x = 0(e 2 / 3 ), together with (|2.41[) and (|2.42 &P for <p in the plateau and transition regions, we estimate 

fL / poo \ / poo \ { pi 



4>dx = £ 2/3 ZT 1/3 Qf $dy^) + £ 2/3 ZT 1/6 Qf ^dy^j ^ ^-jip a dxj +ec J sech 2 (y/2) dy , 

= e 2/3 B- 1/6 J Qdyll + j —S-^fadxj+tec. 
Therefore, the constant c in (|2.41[) satisfies 

" 1/3 , ^--^(T^dy) [! + /(&)] , I(P ) = J^^jiiodx. (2.43) 



c e 



Here B > is defined in (|2.24[> . Since / °° <I> dy > by Theorem and V'o > on x £ (0, f), we get that c < 0. 
Numerically, we compute from (|2.3ip that / °° dy w 1.1857. Alternatively, J(/?o) must be calculated numerically 
from the solution to (|2.42 ap . We remark that the integrand in /(/?o) is well-defined as x — * 0, since although 
u c —>■ 1 as x — > + , we have ipo ~ Pqx/D as x — * + to cancel the apparent singularity in the integrand. 
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Proposition 3 (Nishiura-Ueyama's Condition 2: Dimple eigenfunction) Consider a single-mesa steady-state so- 
lution at the fold point D — D c . Let <f> be the corresponding eigenfunction. For x — 0(e 2 / 3 ), we have 

<t> ~ $ (V/ 6 £- 2 / 3 :r) , $(0) = 1 . 

Here B = (/3q — l)/(2D c ) > and $ (y), defined in \2. 31\) of Theorem^ at s — s c , is a strictly positive function 
that decays at infinity. Alternatively, in an 0(e) region near x = I, we have 

^coe-^sech 2 ^) , 

where cq is the negative constant, independent of e, given in In the outer plateau region < x < I, then 

<f> = 0(e 2 / 3 ) is determined from and this outer approximation for <f> has a unique zero crossing at x = I. 

This establishes the dimple-shape of cj) when e« 1. 



2.3 Universality of the Core Problem 

In this section we show that the core problem can be derived for a class of reaction-diffusion systems that have 
steady-state mesa solutions. On x G [-L.L], we begin by constructing a single mesa steady-state solution for 

u t = e 2 u xx + a(u,v) , av t — Dv xx - v + b(u, v) ; u x (±L, t) = v x (±L, t) = . (2.44) 

We assume that there exists three roots to a(u,v) = on the interval < v < v m at u = 0, u = U-(v), and 
u = m + (w), with < U-(v) < u + (y). Furthermore, we assume that 

a u (0, v) < , a u (u-, v) > , a u (u + , v) < , for < v < v m . (2.45 a) 

We write the two roots u — u±(v) on < v < v m as v = h(u). When v — v m the two roots are assumed to 
coalesce so that u m = w_(u m ) = u + (v m ) and v m — h(u m ). Furthermore, we assume that there exists a unique 
value v c with < v c < v m such that the Maxwell line condition 

lie 

a(u,v c )du — 0, u c = u+(y c ) (2.45 b) 

is satisfied. We also assume that h'(u) < for u > u m and h'(u) > for u < u m . With these assumptions on 
a(u,v), we conclude at the coalescence point that 

°°uu = a U u(u m ,v m ) < 0, a° = a v (u m ,v m ) < 0. (2.45 c) 

For the function b(u, v) in (|2.44|) we will assume that 

6(0, v) = 0; g(u) = h(u) - b[u,h(u)\ < for u > u m . (2.45 rf) 



A specific example of (I2.44[) is the Gierer-Meinhardt model with saturation where a(u, v) = — u+u 2 /[v(l + ku 2 )] 
with k > 0, and b(u,v) — u 2 . For this system we calculate 



u±{v) = 



1 

2kv 



1 ± \J\ - Akv 2 



v = h{u) 



l + ku 2 



h'(u) = -j- 



1 - ki 



(2.46) 



(1 + ku 2 ) 2 

Hence, v m — l/[2yk], u m — 1/Vk, and h'(u) < for u > u m . In Fig. 7(a) we plot u = u± (v) and in Fig. [7(b)1 we 
plot v = h(u). The Maxwell- line condition (|2.45 b\ is satisfied when (cf. |17j ) 



0.4597 

Vk 



u + (v c ) = 



1.515 

Vk 



(2.47) 
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(a) u±(v) (b) v = h(u) 

Figure 7. Left figure: Plot of u±(v) from (|2.46|) for the Gierer-Meinhardt model with saturation parameter k = 
0.25. Right figure: corresponding plot of the inverse function v — h(u) from (|2.46|) 



In addition, we calculate from (12.45 a| that 

9{u) = TTku^ t 1 " u{1 + ku ^ ' (2 ' 48) 

Since g'{u) < for u > u m = \j\fk 1 and g (l/y/jcj < when < k < 4, we have g(u) < for u > u rn when 
< k < 4. 

We now return to the general case under the assumptions (|2.45[1 and we construct a single mesa steady-state 
solution of the type shown Fig. 2] We first derive an expression for the critical value D c of D for which no single 
mesa steady-state solution exists when D < D c . 

Near the interface at x — I we introduce the inner expansion 

u = U (y) + eU 1 (y) + ... , v = V Q + eV 1 (y) + . . . , y = e-\x-l). (2.49) 

From the steady-state problem for (|2.44|) . we obtain 

Ug + a(U o ,V o )=0, Vf=0, (2.50 a) 

U? + a u {U , V )Ui = -a v (U , V )Vi , V{' - . (2.50 6) 



The solution to the leading-order problem is Vo — v c , where v c satisfies (|2.45 &p . and Uo(y) is the unique 
heteroclinic connection satisfying 

f/ (-oo) =u + {v c ) =u c , U {oo)=Q, U {Q)=u c /2. (2.51) 

At next order we obtain that Vi = Vny + V12, for some constants Vn and Vn. The solvability condition for 
(|2.50 b\ determines V\i in terms of Vn as 

/OO /"OO 
a v (U Q ,v c ) U^dy = -Vn / a v (U ,v c ) yU' Q dy . 
-00 J — 00 

Then, by matching to the outer solution for v we obtain Vn — v'(l ). 

The outer problems for v determine v'^). In the mesa region < x < I, where v = h(u), we readily derive 
the following outer problem 

Dv" = g(u), 0<x<l; v(l)=v c , i>'(0)=0. (2.52) 

Here g(u) is defined in (|2.45 The corresponding u is given by u{x) — u + [v(x)] with u(l) — u + (v c ) = u c . We 
require that < v < v m at each x € (0, 1) so that u > u m on x G (0, 1). In contrast, since 6(0, v) — by (|2.45 dl . 
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(2.53) 



Since V\\ is a constant, the solutions to (|2.52p and (|2.53[) are joined by the condition that v'(l~) = v'(l + ). 

The reduction of (12.52)) and (|2.53p to a quadrature relating u(Q) = Uq to the length of the domain L is very 
similar to that done for the Brusselator. We first multiply (|2.52[) by v' — h'(u)u' and integrate to get 



F(u;u ), F(u;u )= / g(s)h'(s)ds 



D- 



Since h'(u) < and g(u) < for u > u m (see (|2.45 d\i ), we obtain that F(u;uq) > for u > uq. By taking the 
negative square root, we calculate 



dv 



du 



^ = -V^ VF(M;Uo)<0 ' d-x=-\D h>(u) >0 - 



(2.54) 



By integrating (|2.54[) with u(0) — u and u(l) = u c = u + (v c ), we obtain a relation between the half-length I of 
the mesa and u$ 

- {Vc) h'(u) J ^2F (u + (v c );u ) , f u +^ g> (t») 



D 



y/2F (u;u ) 



V2F ( M+ ( Uc ); Mo ) 

s[«+(v c )] 



In the outer region I < x < L, we solve (|2 . 53[) to obtain 



v(x) = v c = v x (l + ) 



cosh 


\l- 


x)/Vd 


cosh 


\l- 


i)/Vd 



D 



tanh 



r^2F {u;u )du < 0. (2.55) 



(2.56) 



{L-l)/VD 



By setting v x {l ) — v x (l + ), we obtain that 



v ; = tanh 1 



/ y/2F(u + (v c );uo) \ 



when 



^2F(u + (v c );u ) i 



(2.57) 



Finally, upon combining (|2.55p and (12.57[) . we obtain the following expression relating uq = u(0) to L: 



L , ,_ x ( ^2F{u + (v c )-u a ) \ ^2F (u + (v c );uo) 

-7= - X(uo) = tanh = — 

VD \ v c g[u + {v c )\ 



[g(u)] 



: V / 2F(u;u )du. (2.58) 



Noting that F Uo (u; uq) < for u > uq, a simple calculation shows that x( u o) is a decreasing function of uq when 
uq > u m . Therefore, for the existence of a single mesa steady-state solution, we require that D > D Cl where 
D c ^I?l\ X {u m )\\ ' 

Next, we show that the core problem (|2.26|) determines the local internal layer solution behavior near the origin 
when D = D c . In the this layer near y — we expand 



D = D r 



u = u m + 5u\ + ■••, v = v m + 5 vi H , z = x/8 , 

where 8 <C 1. The nonlinear terms in (12.44[) are calculated as 

a fu 2 \ 

a(u, v)~a° + a° u {u - u m ) + -j L (u- u m ) 2 + a° v (v - v m ) + S 2 ( -^a° uu + a% 



(2.59) 



(2.60) 



Here the superscript denotes the evaluation of partial derivatives of a at u = u m and v = v m . In obtaining (|2.60j) 
we used a — a° = 0. By substituting (|2.59|) and (|2.60p into the steady-state problem for (|2.44[) . and choosing 
<5 = e 2 / 3 , we obtain 

a 

u lzz + -^u\ + a° vVl = , D c v lzz = g(u m ) . (2.61) 
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Here g(u), with g(u m ) < 0, is defined in (|2.45 rip . The solution for v\ is written as 

V X =~A- BZ 2 , 

Then, (|2.6ip for u\ becomes 



uu - 2 -a° v (A + Bz 



0. 



(2.62) 



(2.63) 



From (|2.45 cp we recall that < and a° < 0. Finally, we rescale (|2.63|) by introducing C, /i, and A, by 
u\ = CU, z = fiy, and A — Bfi 2 A. Then, (|2.63|) is transformed precisely to the core problem (|2.26p for U(y) and 
A when 



A* 



1/6 



g-1/6. 



c 



a V a a 

u uu \ u 'uu u 'v 



-1/3 



(2.64) 



By combining these transformations, we obtain the following characterization of the internal layer near the origin: 



u — u r , 



v — v„ 



.2/3 



CU(y) 



2£ 2/3 



"nu \ " , uu u 'v 



-1/3 



s 1/3 c/(y). 



2/3 



X 
7573 



4/3g2/3 
2 



UU V 

1/6 



1/3 



gl/6_ 



(2.65 a) 
(2.65 6) 
(2.65 c) 



Here B is defined in f2"l)2"|) . 

Using the result from Theorem [2] for the core problem (|2.26p , we conclude that the bifurcation diagram near 
the existence threshold of D for a single mesa steady-state solution of (|2.44p has a saddle-node structure. Recall 
that at the saddle-node point U(0) « -0.61512 < and A = A c w -1.46638 < (see Theorem H}. Therefore, 
from (|2.65p . we have u(0) < u m and w(0) — u m > at the saddle-node point, as expected. 



For the Gierer-Meinhardt model with saturation where a(u, v) = — u+u 2 /[v(l + ku 2 )], b(u, v) = 
and v m = l/[2vfc], we calculate that 

a° = -2, al u = -Vk , g(u m ) = 1 - 



2y/k 



Vk. 



with < k < 4 . 



= l/Vk, 



(2.66) 



The existence threshold D c = D c (k) can be computed numerically from (|2.58p for a given domain half-length L. 

Finally, we remark on the local behavior of the time-dependent solution to (|2.44p in the internal layer region. 
If we substitute (|2.59p with u\ = u\{z, t) and v\ — v\{z, t), we readily obtain that 



--2/3,, _,. 1 u uu 2 , 
e u u — ui zz h — — u 1 + a v Vi , 



ae 4/3 v lt = D c v lz 



g(u m ) ■ 



(2.67) 



We then introduce C and r defined by t = e 2 / 3 /j, 2 t and U\ = CU. In this way we obtain, ae 2 fi 2 V\ t 
D c vi zz — g{um)- Thus, v\ is quasi-steady, and D c v\ zz = g{u m ). The corresponding equation for U(d,t) is 



U T = U yy - U 2 + A + y 2 



(2.68) 



If we take A < A c and even initial data U(y, 0) with U (0, 0) < s c , which is below the existence threshold for 
the steady-state core problem, then (|2.68p should exhibit the finite-time blowup U — > —00 as t — > T~ . The local 
structure of the solution near the blowup point y = and r = T is independent of the lower-order terms A + y 2 
in (|2.68p . and is given from [7] as 

,,2 



U(y,T)~-(T- T ) 



-1/2 



1 



1 



y 



(2.69) 



4|log(T-r)| 8(T-r)|log(T-T)|. 
The analysis leading to (|2.68p is, of course, not a valid description of the solution to (|2.44|) when r — * T~ since full 
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nonlinear effects in (|2.44[) must be accounted for near the singularity time. However, this analysis does suggest 
the formation of a large amplitude finger, such as shown in Fig.[lfb), when D is reduced significantly below D c . 



3 The Stability of the Mesa Pattern 

In this section we show that the steady-state if- mesa pattern is stable when D > D c and r — 0. We linearize 
around this steady-state solution by letting 



u (x, t) — u (x) + e xt cf)(x) , w (x, t) = w (x) + e xt ip(x) , 
to obtain (|2.33p . Upon setting r = in (|2.33p , we obtain the eigenvalue problem 

= £ 2 4> xx - f u (u, w)4>- f w (it, w)tp, —<j>= Di/j xx - M , 

a 

where f(u, w) is defined in (|2.7[) . The main result of this section is as follows: 



(3.1) 



Theorem 4 (Nishiura-Ueyama's Condition 3) Consider a symmetric K-mesa steady-state solution as constructed 
in Proposition^ on a domain of length 2KL, with D > D c . When t = 0, the spectrum of \3.1)) admits only real 
eigenvalues with A < 0. 

To show Theorem^] we first reformulate p.l[) as a singular limit eigenvalue problem (SLEP) in the limit e — > 
with a = O (e 2 ) , in order to derive a reduced set of equations, independent of e, for the eigenvalues. The following 
Lemma characterizes this reduced system and its eigenvalues: 



Lemma 5 (SLEP reduction) Let A be an eigenvalue associated with the K-mesa steady-state solution w, u on an 
interval of length 2KL, with an interface at x = I, as described in Proposition^ The leading 2K eigenvalues of 
the eigenvalue problem h3.1\) are of order O (a) with a — 0(e 2 ). These eigenvalues are characterized as follows: 
Define Ai by 

A = aAi , 

where a = 0{£ 2 ), and let u e , u a be the solutions of the differential equation 



£>V"-(/3o + Ai)— ^ = 0, 
w 

satisfying the boundary conditions 

u o (0)=0, u' (l) = l; <(0) = 0, <(/) = ! 

Then Ai satisfies 

--(A) + Ai) = ^, 
where a is one of the eigenvalues of the 2K x 2K matrix 



M 



s 

b 1_ J_ 

S 2d 2d 
J_ _b 1_ 

2d 8 2d 

a 



b 1_ 

" <5 2d 



2d 



J_ 

2d 



S 2d 
a 
8 



(3.2) 



(3.3) 
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The entries of the matrix M are 

a= Me(0 ~"° W , h= ««(0 + «°(0 f S = b*-a\ d = L-l. (3.4) 
The eigenvalues of M are given explicitly by 

i*br Z/ie special case of one mesa, where K = 1, there are two small eigenvalues corresponding to either an even 
or an odd eigenf unction. These eigenvalues satisfy 

{f3 Q + \ l )u e {l) = ^J- and (P + \ 1 ) Uo (1) = ^-1Jl. (3.6) 

This proof of this lemma is given in Appendix A. Here we will use it to prove Theorem |4j 
Proof of Theorem [4j Define /x by 

fx = Ai + Po ■ 

From (|2.10p we obtain on the interval x £ (0, Z) that ^7 = > since u € (l, v2] on this interval. With this 
preliminary result, the proof of Theorem Q] consists of four steps. 

Step 1: Let fi = Ai + 0o and define / (/i) = fj,u a (I) . In this step we will show that the function fi — > uq (I) is 
decreasing whereas / (/i) is increasing for all fi > 0. 

The former claim is easy to show. Indeed, let Ui be a solution of u" — hiU = with m (0) = 0, (I) = 1, 
for i = 1,2, and with < Zii(x) < Zi 2 (x). Then, from the comparison principle, we find that iti > u 2 . Now take 
< [11 < p:%. Applying this comparison principle with hi — an d h% — ^jy^p-, we immediately find that 

u (I; hi) > u (I; fin). 

Next we show the more difficult result that / (/x) is increasing. Define h(x) = ^7 so that uo satisfies 

u" - fih (x) u = , h > ; u a (0) = , < (0 = 1- ( 3 - 7 ) 
Define « and by u = ^ (/ztxo) and w M = gjjU. Then, we readily obtain 

v" — [ihv = fxhuo ; v (0) = , u (Z) = 1 , 
- /xto^ = 2hv ; (0) = , uj, (0 = 0. 

First note that by the maximum principle, u > for all x G (0, I) so that v" — /ihv > in (0, 1) . Now suppose that 
v(l) < 0. Then, by the maximum principle, v < for all x € (0, Z) . But this implies that v'' — fihv fl < inside (0, 1). 
It then follows by the maximum principle that t> M > for all x G (0, Z). Since /' (/i) = v (I) and /" (/x) = (I), 
we conclude that /" (/z) > whenever /' (/x) < 0. It follows that / has no local maximum. Therefore, to complete 
the proof of Step 1, it suffices to show that /' (0) > 0. 

For fj, <C 1, the leading-order solution to ()3.7|1 is u Q (x) ~ x. It follows that / (/z) ~ zzZ as /x — ► so that 
/' (0) = I > 0. This completes the proof of Step 1. 

Step 2: We show that er m i n < a < where 

O-rmn = ~ ( — + \ J , d = L-l. (3.8) 

Vu (Z) a/ 

To show this result, we must establish that < u (Z) < u e (Z), where u Q (Z) = b — a, and u e (I) = b + a from (|3.4p . 
This result follows from a comparison principle, which yields that < u Q (x) < u e (x) for all x € (0, 1]. A simple 
calculation shows that the result a m i n < a < readily follows from (|3.5p upon using < u (Z) < u e (Z). This 
completes the proof of Step 2. 

Step 3: Since a > <7 m i n , we derive that 

- - (ft + Ai) > Q{ji) = t 1 ^— ■ (3.9) 
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u' 

u" D - —[3 u o = ; u (0) = , u' a (I) = 1 . 
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(3.10) 



By differentiating ([2~22]) . we note that u/ satisfies (f37T0|) on [0, 1] with u/ (0) = 0. Therefore, u a (x) = w'(x)/w"(l). 
We then calculate using (|2.22[) that 



•j" (0 = ~(/W2-l 



and w/ (/) 



L-l 



Therefore, for // = /3o, we obtain u Q (I) = (L — Z)/(/3qv2 — 1) and consequently 

P u o (I) L - 



Q((3o) = -r 



when [i = flo . 



Next, by Step 1, we readily find that G(fx) in (|3.9[) is an increasing function of fi whenever /i > 0. Therefore, 



■-(/3q + A 1 )>^(/3 ) = ^ 7 J 
cr V2 



for all fi > flo . 



We conclude that (|3.2[) cannot be satisfied if Ai is real and positive. 

Step 4: To complete the proof of Theorem 01 it suffices to show that all roots Ai to 
do so, we decompose M into the two block-diagonal matrices 



are purely real. To 



M = M 1 +M 2 . 



where 



Mi 



b 


a 


8 


5 


a 


b 


8 


8 



b 


a 


8 


5 


a 


b 


8 


8 



b 


a 


8 


8 


a 


b 


8 


8 



1 


1 


2d 


2d 


1 


1 


2d 


2d 



1 


1 


2d 


2d 


1 


1 


2d 


2d 



We first note that (|3 . 2[) is equivalent to 



Mv = -fi 



L-l 



(3.11) 



where v is an eigenvector corresponding to a and /i = Ai + /3q. Let be the k component of v. Then, upon 
using 5 — b 2 — a 2 , we calculate the inner product as 



Vl a + M a 



-Id- 

"2(6 
C 2 



-!-|V2K-l| +\V2K\ ) + -jj (V1V 2 + V 2 Vi H h + «2K«2i<:-l) 



(3.12 a) 



^— -y - l> 2 | 2 H h |"l>2if-l - f2K| 2 ) - 2 ^ 1 | _ ^ + V 2 \ 2 H h |«2K-1 + «2A'| 2 ) ; 



M Q (/) U e (I) 



(3.12 6) 
(3.12 c) 
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Here and below Ci denotes a non-negative constant that may change from line to line. Similarly, we obtain 

v l M 2 v = -— {\v 2 - H V \v2K-2 - V2K-i\ 2 ^j , (3.13 a) 

= -C 3 . (3.13b) 

Upon premultiplying (|3.11|) by v*, and using (|3.12|) and (|3.13|) . we then divide by fi to obtain 

C\ t C 2 , C-j, 

— — O4 . 



flU e (I) [IU (I) \l 

This equation can be rewritten as 



dfiue (I) + C 2[ iu (0 + C 3 /2 = Ci . (3.14) 

From the expressions for d in (|3.12|) and (|3 . 1 3|1 it follows that at least one of the C\, C%, or C3 is strictly positive 
for any v 7^ 0. Next, we return to the equation for u a , 



11 

Du';-v—u o = 0; u o (0) = 0, <(/) = ! 



We multiply this equation by and integrate the resulting expression by parts to get 



pi pi t 

u^{l) — I \u' \ dx + — / — \u \ dx . 
Jo u Jo w 



Multiplying this expression by /1 we obtain 



liu (I) = [iB 5 + B 6 . 

In a similar way, we derive 

[iu e (I) = JiB 7 + B 8 . 

Here B$, B§, B 7 and B% are strictly positive constants. Substituting these expressions into (|3.14p we conclude 
that 

fl(C 1 B 5 + C 2 B 7 + C 3 )=B, 

Here B is a real constant, and we note that C\B^ + C 2 B 7 + C3 is strictly positive for any 11 ^ 0. Finally, by taking 
the imaginary part of this expression, we get Im (/1) = Im Ai = 0. This concludes the proof of Theorem HI 



4 Discussion 

In [15) the one-dimensional Brusselator model was analyzed with the following scaling, 

TkU t = e k D k u xx + e k A k + u 2 v - (B k +e k )u, v t = e k D k v xx + B k u - u 2 v , 

on the interval x € [0,1] . The assumptions on the parameters were that e k D k <C 1, D k 1, A k = 0(1), 
and B k = O(l). This model is equivalent to (|1.2|) after the change of variables u = au, v = av, t — %t, with 
a = \/B k + ifc, and after dropping the hat notation. The parameters are mapped to 



e k D k \/gfc + gfc n n « 1 

Po = ; D = D k (j ; r 



B k + e k A k r k 

One of the main results in |15j was that a if -mesa configuration with K > 2 is unstable only when 
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Figure 8. Space-time plots showing self-replication of a moving mesa, (a) D = 0.1; (b) D = 0.05. The other 
parameter values are fixed at e = 0.01, (3q — 4, r = 0.1. Critical threshold values for K mesas are Di = 0.65, 
D 2 = 0.16, D 3 = 0.07, £> 4 = 0.04. 

For Dk above this stability threshold a coarsening phenomenon was observed in |15) . This process resulted in 
the annihilation of some mesas over an exponentially long time scale, until eventually the number of mesas was 
decreased sufficiently so that (|4.1j) no longer holds. 

The results in this paper together with [15] provide analytical bounds on D for the existence and stability of 
a steady-state if-mesa pattern. When /3 > V%, (|4.1[) reduces to 



This provides an exponentially large upper bound on D for the stability of K mesas. Roughly speaking, the 
stability of if-mesas when r is sufficiently small is guaranteed on the range 



If D exceeds an exponentially large upper bound, then the number of mesas is diminished through a coarsening 
process. Alternatively, if D is too small, then self-replication is observed until such time that DK 2 is large enough. 

There are several open problems. When the initial data is non-symmetric, it is possible to observe a sequence 
of self-replication events without changing parameter values. An example of this phenomenon is shown in Figure 
[51 In this case, the motion of the mesa itself causes successive replication events until a steady stable state is 
reached. The minimum number of splitting events is given by K = y/ D c /D where D c is given by Proposition [U 
although numerical simulations indicate that the eventual number of mesas is higher than this minimum. Note 
also that a sequence of "abortive" splittings is observed; we speculate that these are connected to the overcrowding 
phenomenon as described in [15] . although the details are unclear. The question of how the motion of the mesa 
affects the splitting behaviour is also open. 

Another open problem is to study the stability and dynamics of equilibrium and quasi-equilibrium mesa patterns 
when r > 0. In [15j . it was shown that when D ^ 1, there is a Hopf bifurcation that occurs for r ~ 1. However, 
the analysis there relied on explicit analytical calculations of the small eigenvalues, which is not possible when 
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D = 0(1). Moreover, it was shown in [15) . under some additional conditions, that a Hopf bifurcation can occur 
leading to a breather-type instability whereby the center of the mesa remains stationary, but its width slowly 
oscillates in time. For D above the self-replication threshold of a steady-state mesa, we suggest that such a 
breather-type instability for r sufficiently large can trigger a dynamic mesa self-replication event if the time- 
oscillating mesa plateau width exceeds its maximum steady-state value. Such a triggering mechanism is explored 
for a reaction-diffusion system with piecewise-linear kinetics in [lOj . In addition, when D = O(l), some numerical 
simulations (not shown) suggest that an oscillatory traveling-wave instability is also possible, whereby the position 
of the center of the mesa oscillates in time, while its width remains constant. 

The second area of open problems is to extend the study of the existence and stabilty of mesa patterns to two or 
higher dimensions. Numerical simulations suggest a slew of possible patterns. One possibility is to study radially 
symmetric patterns in a ball. One can then obtain a blob-like pattern. Self-replication of such a pattern can occur 
as D is decreased sufficiently, and we expect the core problem in two dimensions to be (|2 . 26|) with y replaced by \y\. 
The study of mesa blob-type patterns in an arbitrary two-dimensional domain is also open. Another possibility is 
to trivially extend the one-dimensional mesa pattern into the second dimension. The resulting mesa-stripe pattern 
can exhibit transverse instabilities. This stability problem was examined in [lTj for the Gierer-Meinhardt model 
with saturation in the near-shadow limit where mesa self-replication does not occur. A similar stability analysis 
in the mesa self-replication regime is open. The analysis of these and related problems in two dimensions is the 
subject of future work. 
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Appendix A The SLEP Reduction 

Proof of Lemma [5j We label the interface locations by 

X\- < X\ + < X 2 - < x 2 + < ■ ■ ■ < XK- < x K + , 

as illustrated on Fig. [21(b). We then define / and d by 

I = (x t+ - x^)/2 , d = (L - I) = (.t (j+ i)_ - x i+ )/2 . 

By symmetry I and d are independent of i. 

In the inner region near Xi± , we introduce the inner variables 

4> = $ (y) = $ + e$i + ■ • • , ip = # (y) = *o + e*i + • ■ ■ , A = aAi H , y = e~ l (x- x l± ) , (A.l) 

where a = e 2 ao. Upon substituting (|A.1| into (|3.1|) . we obtain the leading-order system 

*o - A. (0b, W ) $o - f w (U , W ) * = , *o = . 

Here Uq and Wo, satisfying (|2.6p . are given in (|2.9p . We take the + sign for Uq in (|2.Q|) for the inner region near 
x = Xi-, and the — sign in (|2.9p for the region near x — Xi+. By differentiating (|2.6|) with respect to y, we get 

*o = Ci±I/5 (y) , *o = 0, (A.2) 

where Ci±, for i — 1, . . . , K , are constants to be determined. 

Since = 0, we obtain the following problem for <E>i and at next order: 

$1 - fu (U , Wo) $1 - /„ (U , Wo) i = $0 [fun (U , Wq) U x + f uw (U , W ) Wx] , = . 
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In order to match to the outer solution constructed below we require that is a constant. We then multiply the 
equation for <&i by Uq, and integrate the resulting expression by parts, to get 



U' f w {U , Wo) *idy = - $0^0 [fun {U , Wo) U x + f uw (U , W ) W x ] dy . (A.3) 

i J — oo 

To simplify the right-hand side of (|A.3[) , we differentiate the equation for U% in (|2.8p to get 

K - f u (U , Wo) U{ = f w (U , Wo) W{ + U' Q [f uu (U , Wo) Ui + f uw (Z7 , Wo) Wx] = Q, 
Upon multiplying this equation by and integrating the resulting expression by parts, we obtain the identity 

/oo 

f w (Uo,W )Wl$> dy. (A.4) 
- OO 

In (|A.3p . we use (|A.4p , <J> — c i±U' Q , and the facts that and W{ are constants (see (12.121) , to get 

*i / U' U (Uo, W ) dy = c t ±W{ / f w (U , Wo)U' dy . (A.5) 



Since f w = -?7 2 and J^UgU^dy ^ 0, JA3} yields *i = c t ±W{. However, W[ = w'(x t± ) from (l2~T2l . and 
= eip(xi±), where ip(x) is the outer solution for (|3.ip . Therefore, we have the following key relationship: 

ip(%i±) = eci±w'(xi±) . (A. 6) 

Next, we derive an outer equation for ip. In the outer region, defined on the union of the subintervals —KL < 
x < x\- 1 Xi- < x < Xi+ for i = 1, . . . , K, and Xk+ < x < KL, we obtain from (|3.ip the leading-order system 

4> — -h-i> = — jip , Ai0 = Dip xx - j3o(f> ■ 

Ju W 

These equations can be combined to give 

v! 

Dip xx - (Ai + /3 ) — V = . 

w 

To determine the jump condition for ip across x — Xi±, we use use the inner result <p ~ Ci±U^ to derive 

/oo 
U( ) dy = TV2ec l± {X 1 +Po) ■ (A.7) 
-CO 

Therefore, the outer problem for ip on —KL < x < KL is 

Dip" - ([3 + Ai) = -V2e (P + Ai) ( ^ [c 4 +<5 (a: - - c^5 (x - a*.)] J , (A.8) 



with ip'(±KL) = 0. Note that in those outer regions where u — to leading order, we have ^7 = 0. 

Single Mesa: We first analyze (|A.8j) . together with (|A.6[) . for the special case of a single mesa where K = 1. 
For this case xi_ = —I, x\+ = +1, and and x e [— L, L], Then (|A.8|) is equivalent to 



D^"- (A) + Ai)— ^ = 0, i6(-I,I); ^" = 0, ie (l,L)U{-L,-l) ; 



(l + ) - Dip' (r) = -V2e (/3b + Ai) c 1+ , Dip' - Dtp' (-l~) = V2e (j3 + Xi) ci_ 



with ip' (±L) = 0. By solving for ip on x € (/, X) U (— L, — /), this system reduces to 



Dip" — (j3 + Ai) — V> = , a;G(-i,i), (A.9a) 



£>V'(0 = \ / 2e(/3o + Ai)ci + ; (-J) = V2e (/3 + A x ) c x _ . (A.9 6) 
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Figure A 1. Symmetry of ip r , ipi, u a and u e . 

We represent ip in terms of the solutions tpi and tp r to 



with cither 



£>V"-(/3o+Ai)— tp = 0, a: € (-/,/) 



^(-0=0, ^(l) = l or <(-!) = 1, #(0=0. 



(A.lOo) 
(A.lOb) 



We then define a and b by 

a = (-1) , b = i,i{l). 

Since ^ is an even function, we also have 

1p r (—1) = —b , 1p r (0 = ~ a ■ 

In terms of tpi and ip r , the solution for ^ is ip — Mtyl + A r ip r , where ip' (—1) = A r and ip' (+1) = A; 
the boundary conditions for ip in (|A.9j) . we get 



(A.ll a) 

(A.116) 
By satisfying 



4 = ^ (A + Ai) ci+ , 
In terms of this solution, we write the matrix system 

ip t (/) A (i) 1 [ A ; 

tpi (-0 V>r (-0 J [ A r 



A r = ^^(,9o + Ai)ci_ 



V»(0 




. V>(-0 . 





(/?o + Ai 



" 6 


a 




Cl + 


a 


b _ 







To calculate an independent expression for tp(±l) we use the identity (|A.6|) . which states tp(±l) 
To calculate w'(±l), we recall that Dw" = —1 for x £ (—L,—l) and for x £ Q,L). With w'(±L) ■ 
w'(±l) = ±{L - I)) D. Therefore, 



tp(l) 


_(L-l) 


Cl+ 




D 


_c i- 



(A.12) 

C\±EW'{±1). 

0, this gives 



(A.13) 
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b a 




Cl+ 


_(L-Z) 1 


ci+ 


a b 






y/2 (Po + Xi) 


~ c i- 



The eigenvalues of the matrix on the left-hand side of this expression are b ± a. Therefore, Ai must satisfy 

L — I 1 



b±a. 



V2 [fio + Ai) 

Finally, we rewrite (|A.14[) in terms of new functions u (x) and u e {x) defined by 

U = 0; + tp r , U e = Ipl - lp r ■ 

Then, u Q is odd and u e is even, and both satisfy (|A.1Q a[) with the side conditions 

uo(0) = 0, <(0 = 1; <(0) = 0, u' e (l) = l. 
Moreover, by using (|A.11|) and (|A.15|) . we calculate 

u (l) = b — a, u e (I) = b + a . 
Therefore, from (|A.14[) , the eigenvalues must satisfy 



(A.14) 



(A.15) 



(L-l) 1 



u e (I) , or 



(L-l) 1 



U (I) ■ 



V2 (A) + Ai) LW ' ^2 (/?o + Ai) 

The corresponding eigenfunctions are either even or odd. This proves (|3.6[) of Lemma [5] 

General case: We now consider the case of K mesas with K > 1. On each subinterval we solve for tp to obtain 

1p = Aillpu + A ir lp ri , X G [Xi- , Xi+] , 

ip = d + Di (x - x i+ ) , x G [x i+ ,X(i+i)-] U [-.KX, U [x K+ ,KL] , 

where the coefficients An, Aj r , Cj, and A are to be found. The functions ipi, ip r solve I>0" — (Po + Ai) ^7-0 = 
with 

^K-) = o, ^(^i+) = i; ^(xi_) = i, < i (x i+ ) = o. 

Similar to the case of a single mesa, we have 

ipu(xi-) = a, ipu (xi+) = 6, ^ri (a*_) = -6, Vri = -a . 
By satisfying the jump condition for £>?/>' across x = Xi± from (|A.8|) we obtain 

D (0' (x+ ) - V' (xr + )) = A - A !( = - V2e (A> + Ai) c 4+ , (A. 16 a) 

D (0/ (x+_) - V/ (xr_)) = A vr - A-i = V2e (A> + Ai) . (A. 16 6) 

Then, by the continuity of ip across Xi± we get 

aAu - bA„ = Cj_i + 2A-id , bA a - aA ir = , 
where 2d — Xi- — xu-i\+. Then, by using -0' (±KL) = 0, we solve for Cj and A to obtain 



A = D K = , Q = - aAi r , Di = 2d ( a A*+i)/ _ & ^(i+i)r - bA a + aA t 



(A.17) 



Moreover, we calculate 
so that 



ip (xi_) = Ay a - A ir 6, -0 (x i+ ) = A i; & - A ir a . 



An = ^ {bip (xi+) - aip {xi-)) , A ir = ^ (a^ (x i+ ) - &V 0&»-)) , 



(A.18) 
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where 8 = b 2 - a 2 . Therefore, substituting (|A.18|) and (|A.17|) into (|A.16|) . we obtain 



^jj- (A) + Ai) Cj- = |v (^+) - (a*-) 



2(/ 



•\/2e a & 

-p- (Po + Ai) ci_ = ~ip (x 1+ ) - -tp (xi-) 

This system can be written in matrix form as 



V2e{/3 + X^v = Mz, 



2d 

V2e 



2d 

X " £ (A) + Ai) c i+ = ^ - (»»+) - ^V" + ( x (i+i)-) , l<i<K. 



D (A) +Al)CK+ = ~1p{XK-) ~ -^i){x K +) 



Cl- 






-ci+ 




tp(xi + ) 




Z — 








1p (XK-) 









(A.19) 



Here the matrix M is given in (|3.3p . Finally, we use (|A.6[) to calculate ip(xi±) = ±Ci±(L — l)/D. This yields that 
z — e(l — L)v/D. Therefore, (|A.19|) becomes 



Mv= ^L(/3 + Ai)u. 



(A.20) 



The eigenvalue problem (|A.20[) is equivalent to that in Lemma [5] The eigenvalues of M were calculated explicitly 
in |31j , and are given by Q3.5[l . This completes the proof of Lemma [5l 
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